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Abstract 



We study the reconfiguration of the angular resolved photoemission spec- 
trum near M point which occurs in E^S^CaC^Os upon cooling below the 
superconducting transition temperature. Restricting our attention to the case 
of underdoped samples we offer a phenomenological mechanism-independent 
argument which explains the emergence of a peak in the spectrum in terms 
of the normal state pseudogap. All of the basic experimental observations, 
including the peak dispersion, 'dip-and-hump' shape of the superconducting 
state spectrum and appearance of the peak at the temperatures somewhat 
higher than T c , are naturally explained. 



I. INTRODUCTION 

A sharp peak in the angular resolved photoemission spectrum of B^S^CaC^Os 
(BSCCO) is the major subject of this paper This peak emerges from an incoher- 

ent background at temperatures somewhat higher than T c and persists all the way down to 
T = 0. It is located around (tt, 0) (in the units of inverse lattice spacing) in the Brillouin 
zone. Substantial efforts have been invested into uncovering the nature of this peak. It 
was assumed that this might provide us with an important clue about a mechanism of the 
superconducting state. In this paper we will argue that the peak is a consequence of normal 
state phenomenology. Our method will enable us to explain qualitatively the most salient 
features of the peak behavior. 

Unlike overdoped cuprates, it is believed that for the underdoped materials the normal 
state single-particle Green's function can be qualitatively characterized by the presence of 
the pseudogap - a depletion of the spectral density A n k(u) = — 21m G n k(u) near u = for 
k in the vicinity of M point of the Brillouin zone close to the Fermi surface (see fig. lb). The 
origin of the pseudogap is unclear at the moment. Yet, our method does not require such 
knowledge. 

The causality puts some constraints on the Green's function: its real part can be found 
according to the Kramers transformation ||. Then, for any reasonable choice of "pseu- 
dogapped" A n k(cu) the function |ReC7 n k(^)| has two maxima located approximately at the 



edges of the pseudogap (fig. la). We will show that those maxima are responsible for the 
emergence of the 'quasiparticles' in the superconducting state. 

There are several microscopical theories which explain the emergence of the peak below 
T c In the approach presented by J.E. Hirsch charge carriers (holes) are coupled to a 

bosonic bath. If this coupling weakens with the growth of the local hole concentration the 
superconducting pairing driven by the kinetic energy occurs: the kinetic energy benefits from 
the effective reduction of the boson-hole coupling. Another consequence of this boson-hole 
decoupling is the increase of the quasiparticle peak weight. In the paper by M. Eschrig et. 
al. a model of the electrons interacting with a magnetic resonance has been discussed. It is 
claimed that this model reproduces correctly the most salient features of the photoemission 
spectrum. This paper was a development of a previous work by M. Norman, H. Ding 
and collaborators [§,0. The latter is important to us since it has some similarities with 
the present approach. In Ref. [p|,|T0|| the authors concentrated their attention on a step-like 
feature in the electron self-energy which could explain the shape of the spectrum. In order to 
account for such a feature they postulated the existence of a collective mode below T c . This 
mode has to disappear above T c for the model predictions to be consistent with the data. It 
is believed that this mode and the magnetic resonance mentioned above are the same entity. 
Many-body investigation of the phenomenon was undertaken by E. Carlson et. al. They 
assume that in the cuprate materials phase separation takes place. It leads to the formation 
of stripes with every stripe being ID metallic conductor. Using bosonization the authors 
demonstrated that below T c the single-particle spectral density develops a delta-function 
peak at the superconducting gap energy. 

In contrast to these studies, we will deduce this peak without resorting to microscopical 
considerations. It will be shown that the peak in the single-particle spectral density is a 
consequence of the analytical structure of the normal state Green's function. 



II. ANALYTICAL STRUCTURE OF THE SINGLE ELECTRON PROPAGATOR 

In this part of the paper we will calculate the spectral density in the superconducting 
state. Let us derive first an auxiliary relation between the Green's function of an electron and 
a hole. The spin-up electron retarded propag atoris: G c {t) = -z0(t)(c T (t)c{(O) + c{(0)c T (t)). 
The propagator for a spin-up hole: 

G\t) = -ie(t)(c\(t) Cl (o) + Cl (o)c\(t)) = (ie(t)(c\(o) Cl (t) + Cl (t)c\(o))y = (1) 

(im&HQP&SW. + U^c^UJJ- 1 ^)^))* = (-G c (t))* . 

Here U x is the unitary rotation operator. It rotates the spin of the electron by 7r around x 
axis: U~ 1 c a U x = i&af3 c f3 an d U~ l c[JU x = —ia^gct. The ground state is invariant under the 
action of U x . After performing Fourier transformation for the last equality one arrives at: 

G h H = -(G e H)*- (2) 

Assume now that the system at hand is in the normal state but not far from the su- 
perconducting transition. Its retarded Green's function in Nambu representation can be 
written as: 
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r< - ( G nk (uj) \ _ (G nk (uj) > 

GnkM_( o -G ? u(-c)J-l o -G; k (-c)) ^ 



In this form it describes independent propagation of an electron and a hole. Spatial inversion 
symmetry guarantees that G_ k = G k . 

This function G n satisfies Dyson's equation: 

(w - Wo)6nk(w) = i + S nk (a;)G nk (a;). (4) 

Here S n is the retarded self-energy in the normal state: 

Snk(^) = n , , J • ( 5 ) 



o -s; k (-^) 

Next, we change, let's say, doping and drive the system superconducting. New Green's 
function G s satisfies an equation similar to with new self-energy £ sk = £ nk + o~ k where 
a given by: 

Here the relation Exa(^) = (T,2i{—u)))* between two off-diagonal elements of the Nambu 
self-energy has been used. It follows from an equation Gu(lj) = Gl^—u) for anomalous 
propagators which, in turn, can be obtained along the same lines as The point of this 
derivation is to establish the analytical structure of the self-energy in the superconducting 
state. This will be used in our discussion of a sum rule. 

Combining the above equations one gets for G s the expression: 

G sk =(G^-a k y\ (7) 



III. GREEN'S FUNCTION IN THE SUPERCONDUCTING STATE 

The equations derived in the previous section are completely general. In order to proceed 
further we need to make a simplifying assumption: near the Fermi surface electronic {u < 0) 
and hole (u > 0) parts of the spectral function are symmetric for small |a>|. One can express 
that as: 

A k (u) = A k (-u) (8) 

for |co>| < uj g where u g is the size of the pseudogap. Exact particle-hole symmetry is absent 
in the cuprates. We would like to argue, however, that (|) holds at least approximately for 
k close to the Fermi surface. Firstly, the tunneling data ]TlJ] shows that the pseudogap is 
symmetric with respect to voltage polarity change. One might think that this tells little 
about validity of (H) since the tunneling differential conductance measures k-integrated 
spectral density. However, Ding et al. (inset of fig. lb of Ref. 0) demonstrated that scanning 
tunneling microscope (STM) spectrum multiplied by the Fermi distribution function at 
appropriate temperature coincides with the angular- resolved spectral density at k = (ir, 0). 
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This means that STM spectrum equal to A k (u) for uj < 0. It is natural to expect that the 
spectral density coincide with the STM conductivity for uo > as well. Thus, the symmetry 
of the STM data suggests the symmetry of the photoemission spectrum. Although, it is 
not our goal to construct a theory of tunneling into cuprates we may speculate that this 
agreement between the tunneling conductance and the photoemission spectrum is not a 
coincidence but a consequence of the inter-plane hopping matrix element k-dependence: as 
discussed before [12]-^] the inter-plane hopping matrix element is at its maximum for the 



electrons with a momentum parallel to (ir, 0) and it is vanishes along Brillouin zone diagonal. 
Therefore, it is conceivable that the most of the tunneling current is carried by the electrons 
near (vr,0). Secondly, we would like to direct our attention toward the discussion of Ref. 
|15|1 . In the latter paper the angular resolved spectrum at k = kp was analyzed with the 



help of a momentum distribution sum rule. Theoretically, it was proved that an expression 
1 — / °° duj tanh(cj / 2T) (A kF (u) — A kF (— u)) has to be temperature independent provided that 
(§) is true for small \u\ |TB[ . This expression can be easily extracted from the data since 
it is proportional to the integral over u of the photoemission spectrum. The spectrum for 
k near (71*, 0) was integrated, and the integral showed no dependence on the temperature, 
supporting the assumption of the symmetry. Finally, we would like to mention two other 
papers ]TT| , [TE| ] where (^) was assumed to be true. Thus, we believe we have enough evidence 



to treat (|8|) as a reasonable approximation for small \u\. The structure of the background 
for |co>| > uj g is of little interest to us since the positions of |Re G n \ maxima are insensitive 
to the background. For the purpose of simplicity we will assume that the background is 
also symmetric. That is, (||) is true for any u. In such a case Green's function satisfies 
G k (-u) = -G* k (uj). 

A consequence of the symmetry is a constraint on the self-energy in the superconducting 
phase: 

where is a real number and the functions fi and a A satisfy the following: 

^(-cu) = and a A (-cu) = (a A (cu))*. (10) 

In terms of these functions the superconducting spectral density equals to: 

A sk = -21m { ^- (GnOV 1 . (n) 

Let us simplify this equation. If the system is close to the transition point then a is small 
and the denominator of (|TT|) is close to unity everywhere except the edges of the pseudogap. 
At those points as one can see from fig. la (Re G n ) 2 is big and one should be more careful. 
If we are away from the Brillouin zone diagonal where a A is zero then to the second order 
in cr A and to the first order in \i = O^a^l 2 ) we find that the denominator is equal to 
1 — (a A ) 2 (G n ) 2 . In this expression the terms like G n \i have been dropped. Although, they 
are of the order of \cr A \ 2 we presume that they are smaller than (Re G n ) 2 \a A \ 2 since the 
latter is proportional to the second power of Re G n which is big at the pseudogap edges. 
This gives for the spectral function: 
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* r / nk 2 ( ^ 5 , (12) 

where we omit (G n ) 2 n term in the numerator assuming that it is small compare to G n . 
The last formula is the key to understanding the appearance of the 'quasiparticles' below 
T c . At the edges of the pseudogap (Re G n ) 2 peaks up, the denominator gets smaller and 

the speetra! taction delps a ship rJJL. This * true provided thft Re{(.-f} is 

positive at the edges of the pseudogap |ID . 

IV. SPECTRUM PROPERTIES 

Now we would like to analyze (|ll|) and (|i~2| ) to obtain simple properties of the 'quasipar- 
ticles' which can be compared against the experiment. 

First, it is clear that the position of the peak in the frequency domain is determined by 
the pseudogap edge. Therefore, we expect to see strong correlation between the size of the 
pseudogap and the frequency (binding energy) position of the peak. This is precisely the 
effect reported in |TJ where temperature dependence of the peak frequency was compared 
against that of the pseudogap size. Fig.2C of that paper shows that the peak position 
traces the size of the pseudogap. Similar phenomena is seen for the peak dispersion (binding 
energy vs. k). As reported in |2tJ the pseudogap at (tc, 0) is virtually k-independent. This 



is consistent with the weak dispersion of the peak itself |2jJ . 

In order to obtain the shape of the spectral function at k = (71", 0) we performed simple 
numerical study. We modeled /x, a A and A n as 

H = 

R ea A = ar£exp(-(u/u A ) 2 ) (13) 

Aik = A (arctan ((a; — u> 9 ) /uj e ) — arctan {(u + u g ) /uj e ) + tc) exp (— |u;/W|) . 

The energy u e determines how abrupt the edges of the pseudogap, W mimics the bandwidth. 
The constant A has to be determined from the sum rule ||: 

L ^ Mu) = L (14) 

The parameter u A describes the frequency dependence of the anomalous self-energy. It is 
commonly put equal to infinity. This makes the anomalous self-energy frequency indepen- 
dent. Below we will discuss both cases of finite and infinite uj a . 

The imaginary part of cr" 4 and the real part of the Green's function can be found with 
the help of numerical Kramers transformation. Then the spectral function in the super- 



conducting state is determined using (|TT|). The resultant spectra for different parameter 
values are presented on fig. 2 and fig. 3. These two figures correspond to different values of 
uj a {uj A /ujg = 15 for fig.2 and uj A /ujg = 7 for fig. 3). On each figure the top (bottom) row 
corresponds to a A juj g = 25 (a A /uj g = 15). The value of u e /u g in the left (right) column is 
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0.3 (0.1). We see clear 'dip-and-hump' structure of the spectrum for big values of <j a . When 
this quantity is small we find only a tiny bump or even a kink instead of a well-developed 
peak. Such situation seems to be realized in Pb-doped BSCO compound P2"| . This angular 
resolved photoemission study revealed no peak in the spectrum (fig. 3, right panel of the 
latter reference). Authors claim that their energy resolution is high enough to observe the 
peak should it be present in the spectrum. 

By comparing different spectra from fig. 2 and 3 we can discuss the dependence of the 
spectrum shape on the model parameters. Two parameters which affect the spectrum the 
most are a A /u) g and uj e /uj g . The effect of a A is obvious: the bigger the anomalous self- 
energy the higher the peak. The ratio u) e /u g regulates the height of the peaks of Re G n . 
This height is proportional to — \n{uo e / uj g ). Thus, the smaller this ratio the higher the peak in 
the superconducting spectrum. The value of u A specifies roughly the position of the hump. 
For u A = oo the hump is absent from the spectrum - it is moved to infinite frequency. 

Another interesting property of our superconducting state spectrum is its compliance 
with the sum rule (|14|). Mathematical proof of this fact is given in Appendix. Visually, one 
can notice from figures 2 and 3 that the depletion of the spectral weight in the dip is roughly 
compensated by its enhancement at the hump. Experimentally, however, it is seen that the 
weight from the dip goes into the peak. The reason for this discrepancy is the lack of exact 
knowledge about a. Particular, we are unable to identify any first principal based constraint 
which would provide us with any information about [i. To keep our consideration as simple 
as possible we chose to put fi equal to zero. 

As we stressed the equality (§) is only an approximation. Thus, an interesting question 
worth discussing is how deviations from (|8|) affect the peak in the superconducting state. 
To address this issue we model the electron spectral function as 

A c nk (uj)=A n]i (uj + Auj), (15) 

where A n k(cu) is given by (fL3|). In this case the pseudogap for electrons is bigger then the 
pseudogap for holes by the amount of 2Au. The spectral function can be found according 
to 

A e k = —2lm I ^ nk I (lg) 

I 1 _ (^k ) Gnk^nk J 

The resulting function for Atu = 0.2u g and <j a (uj) = 25u g is plotted on fig.4. As one can 
see the only effect of the asymmetry is uneven hight of the peaks for positive and negative 
frequencies. The reason for this becomes clear if one notice that (|16|) may be approximated 
in the manner similar to (|T2|): 

A e A ~ 2 A " k f T . (17) 

1-KfRe {G* k G-k} 

In order to demonstrate the quality of this approximation we plotted this expression on 
fig.4. The denominator in ( fT^) is symmetric with respect to the frequency sign change. It 
has two minima located approximately at ±uj g . Since A^ (u g ) > A e n {— uj g ) the peak at uj < 
is smaller. It is now clear that the superconducting state spectrum is rather robust with 
respect to small deviations from (§|). 
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To understand why the peak appears at the temperatures higher than T c it is enough to 
imagine that in the cuprates there is a temperature window Ti ocal > T > T c where the phase 
coherence exists over some finite time scale r(T), r — > oo as T — > T c . Experimental evidence 
of this was produced by Speaking more technically, in this temperature window one 
can use (Jfty with non-zero a A and fixed (f) to describe the single-particle dynamics on the 
time scale less than r. For time periods bigger than r the dynamics of the phase <ft cannot 
be neglected anymore. If r is big enough one can simply average the propagator over the 
phase fluctuations. For infinitely long time the anomalous propagator always vanishes above 
T c due to the phase dynamics. However, the peak in the photoemission spectrum does not 
vanish since, as it follows from flTT|), it is independent of cf). 



V. DISCUSSION 

We study the 'quasiparticle' peak in the superconducting state of the underdoped 
BSCCO. The main conclusion of the present work is that the peak can be explained on the 
basis of the normal state phenomenology. Our analysis is based on three premises. First, we 
assume that the single particle spectrum has a pseudogap (fig. lb). The most important con- 
sequence of that is the peaks of the function |Re G n (uS) \ at the edges of the pseudogap. Next, 
we assume that there is approximate symmetry between positive-frequency and negative- 
frequency parts of the spectrum (||). This guarantees that the function —Re {G n (u)G n (— u)} 
has only one maximum at u < 0. We also shown that some violation of ([8]) is not fatal for 
the qualitative structure of the superconducting state spectrum. Finally, there is the third 
assumption we made implicitly: the single-particle anomalous self-energy a A (u) does not 
go to zero at the edges of the pseudogap. If a A does go to zero then the peak in the 
superconducting phase could be completely suppressed (see (|12D). 

In order to explain the emergence of the peak at T > T c we accepted another postulate. 
It was agreed that (||) local order parameter is established and it is permissible to use @ 
with fluctuating to describe the single-particle dynamics above T c . 

The fact that the peak can be derived from the normal state properties suggests that it 
is not a fundamental object. Moreover, depending on the circumstances it may be absent 
or replaced by a less prominent feature p2| . 



It is also can be surmised that this phenomena may not be limited to the cuprate super- 
conductors only. The peak could be present in the superconducting state of any material 
which has a pseudogap in the normal state. 

We did not try to fit any experimental data by our curves from fig.2 and fig.3. It is 
not possible owing to the lack of knowledge about exact form of functions flT3|). We also 
have no independent information about the values of ui e , u A , a^. Another argument against 
feasibility of the spectrum fit is a problem with the experimental data themselves. It is known 
that the photoemission spectrum is a convolution of the single-electron Green's function and 
a photoemission matrix element. It has been assumed that this matrix element is a constant, 
at least in a relevant range of parameters. Thus, relative intensities of different features of 
the spectrum were believed to be free of any influence extrinsic to the electron dynamics. 
However, it was shown recently that the matrix element has very non-trivial dependence on 
the photon energy and the binding energy of the photoelectrons p^ j. In such a situation 
the significance of an intensity fit is substantially devaluated. 



7 



To conclude, we study the peak in the photoemission spectrum of BSCCO. It was shown 
that this peak is a consequence of the normal state pseudogap. We derived the most basic 
properties of the peak and demonstrated that they are in qualitative agreement with the 
experimental data available. 
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VII. APPENDIX 



In this Appendix we will prove that under rather mild assumptions on G n and cr A the 
spectral function ([11]) satisfies (|14|). In our derivation we will assume that fi = 0. This 



limitation is not crucial. 

An essential part of our proof is the analyticity of both G n {uj) and a^iu) for complex 
values of u, Im u > (upper half-plane). First, we establish asymptotical behavior of G n at 
large \u>\. If |u;| 3> W the details of the spectral function A(u) are irrelevant and it can be 
viewed as a delta function times 2ir. This means that at large \u>\, Imw > 

GnM = / T, 7 - -• 18 

J-oo 2tt u — U) uo 



We see, that due to the sum rule fll4|) the residue of G n at infinity is exactly unity. Next, 
for any finite uj a the same reasoning can be applied to a A with the result: 

a A {u) = CKM" 1 ) (19) 

for large |u;|. We do not try to determine the residue of a A at intinity since it is of no interst 
to us. The asymptotical behavior combined with the absence of any singularity in the upper 
half-plane of complex uj implies that both G n and a A are bounded: 

\a A {u)\<C ai (20) 
\G n {u)\<C G . (21) 

where C a £ are some real positive constants. Now we are in position to prove that: 

/+oo 

du (G s {u) - G n (u)) = 0, (22) 
-oo 

Gs{u) = i-(4))Vh) 2 ' (23) 

provided that the product CQC a < 1. This inequality guarantees that the denominator of 
G s {uS) never goes to zero. Therefore, the integrand of ( p2[ ) has no singularities for Im uj > 0. 
This means that we can transform the integration contour into a semi-circle in the upper 
half-plane. The integral over the semi-circle of infinite radius equals to the residue at infinity 
times (— 7r). We already know that this residue is unity for G n . It is also unity for G s : at 
large |a>| the denominator is non-singular and approaches 1, thus, the leading asymptotic 
behavior of G s is the same as that of G n . The residue of the integrand is zero which proves 
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FIG. 1. Qualitative structure of the normal state Green's function for different values of io e /uj g . 
The real part of G n {uj) (panel (a)) can be obtained by applying Kramers transformation to the 
imaginary part, panel (b). The solid line (dash-line) corresponds to LO e /uj g = 0.1 {uj e /uj g = 0.3). 





FIG. 2. Superconducting state spectral function (A s (u) vs. co/uJ g ) for W/co g = 50, u /co g = 15 
and different values of cr^/uj g and u e /oj g . The left (right) column corresponds to oJ e /u) g = 0.3 
{uj e /uj g = 0.1). The top (bottom) row corresponds to a^/uj g = 25 (a~^/uj g = 15). The normal state 
spectrum is given in a dash-line. The spectral functions are multiplied by the Fermi distribution 
function with T = 0.2uj g . 
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FIG. 3. Superconducting state spectral function (A,(u;) us. u)/u g ) for W/ui g = 50, uj^/ujg = 7 
and different values of a^/u> g and u! e /ui g . The left (right) column corresponds to u! e /u> g = 0.3 
(u) e /u>g = 0.1). The top (bottom) row corresponds to a A juj g = 25 {a^/ujg = 15). The normal state 
spectrum is given in a dash-line. The spectral functions are multiplied by the Fermi distribution 
function with T = 0.2cj„. 
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FIG. 4. The effect of weak particle-hole assymetry on the superconducting spectral function. 
Normal state spectrum is shifted to the left by amount of Alo = 0.2to g . The values of parameters 
are: W/oj 3 = 50, oj" 4 - = do, oJ e /u g = 0.1, a^/ujg = 25 The normal state spectrum is given in a 
dash-line, the approximation (|l~7|) is plotted in a dotted line. 
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